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Abstract 

Two novel classes of many-body models with nonlinear interac- 
tions "of goldfish type" are introduced. They are solvable provided 
the initial data satisfy a single constraint (in one case; in the other, 
two constraints): i. e., for such initial data the solution of their initial- 
value problem can be achieved via algebraic operations, such as finding 
the eigenvalues of given matrices or equivalently the zeros of known 
polynomials. Entirely isochronous versions of some of these models 
are also exhibited: i.e., versions of these models whose nonsingular 
solutions are all completely periodic with the same period. 
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1 Introduction and main results 



Long time ago the possibility was noticed of using the nonlinear mapping 
among the zeros and the coefficients of a polynomial in order to identify 
solvable (classes of) many-body problems, characterized by nonlinear equa- 
tions of motions of Newtonian type (with one-body and two-body, generally 
velocity-dependent, forces) [1,2]. [Terminology: here and hereafter we de- 
note as solvable any problem whose solution can be achieved by algebraic 
operations, such as finding the zeros of polynomials]. The starting point of 
these developments is a linear PDE such as 

tptt ~ K - (N - 1) a 6 z] ip t + [a 2 + a 3 z-2 (N - 1) a 10 z 2 ] ip z 
- (a 4 + a 5 z + a 6 z 2 ) i/j zt + (a 7 + a$z + a g z 2 + aio-2 3 ) ^ zz 

-N [a 3 + (N-l) (a 9 -a 10 z)} i/; = . (1) 

[Notation: hereafter N is an arbitrary positive integer (N > 2), the symbols 
a,j and (see below) Aj,Bj denote (a priori arbitrary) constants, and sub- 
scripted variables denote partial differentiations of the dependent variable 
ip(z,t); note incidentally that this PDE coincides, up to trivial notational 
changes, with eq. (2.3.3-1) of Ref. [2]]. One then notes that this PDE ad- 
mits solutions which are (for all time) a monic polynomial of degree N in 
z, and introduces the N zeros z n (t) and the iV coefficients c m {t) of such a 
polynomial solution, 

N 

71=1 

It is then easily seen that the fact that ip (z, t) satisfies the linear PDE (JTJ) 
implies that the coefficients c m (t) evolve according to the following system 
of linear ODEs [2]: 

c m - (N + 1 - m) a 4 c m _i - [(N - m) a 5 + aj c m + m a 6 c m+1 

+ (N + l-m) (N + 2 - m) a 7 c m _ 2 + (N + 1 - m) [(N - m) a 8 + a 2 ] c m _i 

—m [(2 N — 1 — m) a 9 + a 3 ] c m + m (m + 1) ai c m+ i = . (3) 



= ^ + ^ Cm (t)^— . (2) 

m=l 
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Likewise, it can be shown (see Ref. [2], or Section 2 below) that the zeros 
z n (t) evolve according to the following system of nonlinear PDEs 

z n = a 1 z n + a 2 + a 3 z n - 2 (N - 1) a 10 

N 

4" ^ ^ •[ {z n £ m ) [2 Z n Z m + {z n 4~ Zffi) (0,4 + G5 Zn) 
m=l 

+a 6 (i„ z m + i m z n ) z n + 2 (a 7 + a 8 2 n + a 9 2^ + ai ^)] } . (4) 

Notation and terminology: hereafter indices such as m, n range from 1 to N 
(unless otherwise indicated), superimposed dots denote differentiations with 
respect to the time t, and the N zeros z n (t) are interpreted as the coordi- 
nates of iV "point particles" evolving according to the "Newtonian equations 
of motion" (J3|) featuring one-body and two-body velocity dependent forces; 
and an analogous interpretation shall be given in all analogous cases below. 
Unless otherwise explicitly specified, we hereafter assume that all quantities 
under consideration (except the time t) are generally complex numbers, so 
that the motions of the particles characterized by the coordinates z n (t) take 
place in the complex z-plane; of course such motions can easily be reinter- 
preted as taking place in a "more physical" real plane, but we shall not 
devote space to discuss this well-known aspect (see, for instance, Chapter 4 
of Ref. [2], entitled "Solvable and or integrable many-body problems in the 
plane, obtained by complexification"). 

Clearly, since the parameters aj are (arbitrary) constants, the solution 
of the system of linear ODEs ([3]) satisfied by the coefficients c m (t) can be 
achieved by purely algebraic operations; and once the N coefficients c m {t) of 
the polynomial ip{ z yt) have been obtained, the computation of the iV zeros 
z n {t) of this polynomial, see ([2]), is also a purely algebraic task. Hence the 
determination of the N coordinates z n (t) whose time evolution is character- 
ized by the Newtonian equations of motion (jlj) can be achieved via purely 
algebraic operations: this iV-body problem is solvable. 

The 10-parameter class of solvable iV-body models (T4J) is vast, and it in- 
cludes several interesting cases [2] . Even the very simplest case, characterized 
by the vanishing of all the parameters aj, is quite remarkable: in this case 
the Newtonian equations of motion read simply 

N . 

z n — — 1 i° a J 

m=l Zn Zm 
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and the solution of the initial-value problem for this problem is given by the 
following simple proposition: the N coordinates z n (t) are the N zeros of the 
following algebraic equation in z, 



Note that, by multiplying this equation by the product ]~[ [z — z n (0)] , it 



becomes a polynomial equation, with time-dependent coefficients, of degree 
TV in z, which indeed generally has N solutions, at least if one allows the 
coordinates z n to be complex numbers, hence their time-evolution to take 
place in the complex z-plane; indeed such an evolution can also be reinter- 
preted as the motion of N point particles on a real ( "physical" ) plane, and its 
phenomenology is quite amusing, see Section 2.4.2 of Ref. [2], entitled "The 
simplest model: explicit solution (the game of musical chairs), Hamiltonian 
structure" . Because of these remarkable features, as well as the neat form of 
the Newtonian equations of motion fl5a|) (which also happen to be Hamilto- 
nian), this iV-body model was given the honorary title of "goldfish" [3], and 
this attribute was then extended to several solvable models featuring in the 
right-hand side of their Newtonian equations of motion ( "acceleration equal 
force") terms such as that appearing in the right-hand side of (I5al) : note that 
these solvable models "of goldfish type" include now the models belonging 
to the 10-parameter class (jlj) and several models outside of this class (for an 
updated review of such models see [4]). 

In the present paper we exhibit two classes of solvable iV-body models 
"of goldfish type" which are, to the best of our knowledge, new, hence, in 
particular, not included in the class (jl]). These novel models are however 
solvable only provided the particle coordinates z n (t) satisfy additional con- 
ditions, which must of course be compatible with the time evolution, so that 
it is sufficient that they be satisfied, in the context of the initial-value prob- 
lem, by the initial data. This limitation could be considered as a drawback of 
these new models, or instead as a feature that adds to their interest in the 
context of mathematical physics. The hunch that such models should exist 
was suggested to us by recent results [7-11] concerning certain exceptional 
polynomial subspaces, as tersely reviewed in Appendix A. 

The first class of these models is characterized by the following (6-parameter) 




(5b) 



N 



n=l 
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Newtonian equations of motion: 

2 Ad 



A 1 z n + A 2 + A 3 

N 



-\- ^ ^ {%"n %ni) [2 Z n Z m -\- A^ (yZ n -\- 2 m ) Z n 
m=l 

+2 (A*- A 2 z n + A 6 z n 2 )]} . (6) 

Here the novel element - not included in the class (jl]) — is that associated 
with the constant A4. In the following section we show that these equations 
of motion entail that the coefficients c m (t) - related to the coordinates z n {t) 
via the mapping (j2J) - evolve according to the following system of linear 
ODEs: 

c m - [{N - m) A 5 + A x \ c m + (N - 1 - m) (N + 2 - m) A 4 c m _ 2 
-(N + l-m) (N - 1 - m) A 2 c m _ x 

-m [(2JV-l-m) A 6 + A 3 ] c m = , (7) 

of course with Co = 1 and c m = for m > N and m < (consistently with (|2J), 
and also with this system of ODEs). Note that the linear character of this 
system of ODEs entails the solvable character of the Newtonian equations of 
motion (jSj). But these findings are only valid if the coordinates z n (t) satisfy 
the constraint 

and correspondingly (see below) the coefficients c m {t) satisfy the constraint 

c*-i(t) = . (9) 

The treatment of Section 2 implies that these constraints, (jSj) respectively 
are indeed compatible with the evolution equations (J5J) respectively (j7|); 
hence they are automatically satisfied, in the context of the initial-value 
problem, provided the initial data for the iV-body problem (jSj) satisfy the 
following two conditions: 

N 1 N ■ (n\ 

y^r = o, y^^ = o, (10) 
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and correspondingly the initial data for the linear system of ODEs ([7j) satisfy 
the two conditions 

CN-M = c N ^(0) =0 . (11) 

Before turning to the second class of new models, let us display an equiv- 
alent avatar of the Newtonian equations of motion 

N — 2 \ 1 A 
A 1 + ^—A 5 ) z n -(N-2) A 2 + [A 3 + (N-2) A 6 ] z n 4 - 



A N ( A 

H — ^- Z + Aq Z + < (z n — z m ) \2 z n z m H — — (i n + z m ) (z n + z m ) 

m=l ^ 

+2 A 4 - A 2 {z n + z m ) + A 6 (zl + z m )]} , (12) 
where 

N 

Z(f) . (13) 

n=l 

[This second version, (Tl2l with (TT3l) . has the merit to immediately yield (by 
summing over n from 1 to N - using and the vanishing of the double sum 
in the right-hand side due to the antisymmetry of the summand under the 
exchange of the dummy indices m and n) the linear ODE 

Z=[A 1 + (N- 1) As] Z - N (N - 2) A 2 + [A 3 + 2 (JV - 1) A 6 ] Z ; 

which is consistent, via the relation 

ci(t) = ~Z(t) (14) 

implied by (j2J), with the equation satisfied by c\{t) that obtains by setting 
m — 1 in ([7]). On the other hand the first version, ([?)]) , has the merit to 
feature only one-body and two-body forces]. 

The second class of new models is characterized by the following (7-parameter) 
Newtonian equations of motion: 

z n = B lZn -(N-l) B 2 z n -2 (N-l) B 3 z* + B 4 ^ 

Zn 

N 

+ { ( z n — Zm) 1 [2 Z n Z m + (z n + Z m ) (_B 4 + B 5 Z n ) 

m=l 

+ B 6 (z n z m + z m z n ) z n + 2 (B 7 z n + B 2 zl + B 3 z^)] } . (15) 

6 



Clearly the novel element (in (1151) relative to (jlj)) is now that associated 
with the constant B^. In the following section we show that these equations 
of motion, (fT5l) . entail that the coefficients c m (t) - related to the coordinates 
z n (t) via the mapping (jSj) - evolve as follows: 

c m -(N - m) B A c m _x - [(N -m) B 5 + B x ] c m + mB 6 c m+1 

+ (N + 1 - m) (N -m) B 7 c m _i 
— m (N — m) B 2 c m + m (m+ 1) B 3 c m+ i = , (16) 

entailing the solvable character of the Newtonian equations of motion (jH]). 
But this conclusion is valid only if the coordinates z n (t) satisfy the constraint 

and correspondingly (see below) the coefficients c m {t) satisfy the constraint 

c N {t) = . (18) 

The treatment of Section 2 implies that these constraints are indeed compat- 
ible with the evolution equations (115j) and (TIB]) ; hence they are automatically 
satisfied, in the context of the initial-value problem, provided the initial data 
for the iV-body problem f|T5|) satisfy the following single condition: 

£St-°- (19) 

and correspondingly the initial data for the linear system of ODEs ([TBI) 
satisfy the single condition 

c N {0) = . (20) 

The paper is organized as follows. In Section 2 the results reported above 
are proven. In Sections 3 and 4 we discuss tersely the behavior of these two 
solvable iV-body problems, and in Section 5 entirely isochronous versions of 
these novel models are identified: they are characterized by the fact that all 
their nonsingular solutions are completely periodic (namely, periodic in all 
their dependent variables z n (t)) with the same basic period (or possibly with 
a - generally not too large [12] - integer multiple of it). The last Section 6 
outlines tersely further developments, to be reported in future publications. 
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The paper is completed by two Appendices: in the first one certain relevant 
results concerning exceptional polynomial subspaces are tersely reviewed, in 
the second one certain computations are confined whose treatment in the 
body of the paper would disrupt the flow of the presentation. 



2 Proofs 



The procedure to arrive at the results reported above is by now textbook 
material [2,4], hence our treatment here can be terse; although the new twist 
should be emphasized, implying that the solvable character of the new models 
applies only provided their time evolution is somewhat restricted - as already 
indicated in the preceding section. The starting point of our treatment is the 
linear PDE (but see the Remark below) 

iptt ~ K -{N — l)z a 6 ] ip t + [a 2 + a 3 z - 2 (N — 1) a 10 ) ip 2 



(a 4 + a 5 z + a 6 z 2 ) ip zt + (a 7 + a 8 z + a 9 z 2 + a w 
N [a 3 + (N - 1) (a 9 - a 10 z)\ i/> 



5 ) 1>, 



flu 

z 



, c N 

Ipt 

cn 



aw 

z 



cn 



(21) 



[Notation: the symbols aj denote again (a priori arbitrary) constants, sub- 
scripted variables denote partial differentiations, and the dependent variable 
ip(z, t), as well as the coefficients cjv(i) and C/v-i(t), are characterized by (T5])]. 

Remark. The consistency of this evolution PDE, (12TT) . with the fact that the 
dependent variable ip{z, t) is the monic polynomial (j2J) of degree N m z shall 
be clear from what follows. Note however that this evolution equation, 02ip . 
is a linear PDE iff the two constants an and a\2 vanish (in which case it 
coincides with (JT])); otherwise it is in fact a nonlinear functional equation, 
as implied by the formulas (clearly entailed by ([2])) 

c N (t) =V(0,t) , cjv-i(t) = V 2 (0,t) • (22) 
Via (El), this evolution equation (|2TT) entails^ the following (systems of) 



1 The derivation of the first of these two systems of ODEs is a rather trivial exercise; 
the derivation of the second system is as well trivial, but readers who find this second task 
exceedingly painstaking are advised to use the formulas given in Section 2.3.2 of Rcf. [2], 
or (even more conveniently) in Appendix A of Ref. [4] . 
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evolution ODEs for the JV coefficients c m (t) and for the JV zeros z n (t): 

c m - (N + 1 - m) a 4 c m -i - [(JV -to) a 5 + Oi] c m + m a 6 c m+i 
+ (AT + 1 - m ) (N + 2 - to) a 7 c m _ 2 + (JV + 1 - to) [(JV - to) a 8 + a 2 ] c m _x 

-to [(2 JV - 1 - to) a 9 + a 3 ] c m + to (to + 1) a w c m +i 

1 ,(23) 



Cat 

«11 I C m _i C m -1 | + Ol2 

Cat 



(JV + 2 - to) c m _ 2 - <^-i. 



Cat 



z n = ai 2„ + a 2 + a 3 z n - 2 (JV - f ) a w z n + a u 

N 

+ { (^n — z m) 1 [2 i n z m + (i n + i m ) (a 4 + 05 Z n ) 

m=l 
m^n 

+ a 6 (i n z m + z m z n ) z n + 2 (a 7 + a 8 z„ + «9 ^ + aio z n )] } .(24) 

The observation that the system (123|) is consistent with Co = 1 and c m = 
for m > JV and m < confirms the consistency with the evolution equation 
(|2"T|) of the assumption that 4>{ z it) be a monic polynomial of degree JV in z, 
see ©. 

In order to guarantee that these Newtonian equations of motion, (T24|) . 
be solvable, we must make sure that the corresponding evolution equations 
( f23|) be effectively linear. A straightforward way to achieve this goal is to 
set to zero the two parameters an and a i2 ; but this simply reproduces the 
previously known (class of) solvable JV-body models (HI). There are however 
two other possibilities, which do yield two new classes of integrable systems. 

The first possibility is to set an = (but ai 2 ^ 0) and to then require 
that the system of evolution equations (|23|) be consistent with the constraint 
(J9j). It is indeed plain that the system (1231) becomes then linear: 

c m - (JV + 1 - to) a 4 Cm-i - [(JV - to) a 5 + ai] c m + m a 6 c m +i 
+ (JV + 1 - to) (JV + 2 - to) a 7 c m _ 2 + (JV + 1 - to) [(JV - to) a 8 + a 2 ] c m _x 

-m [(2 JV - 1 — to) a 9 + a 3 ] c m + m (to + 1) a w c m +i 

= a 12 (JV + 2-m) c m _ 2 ,(25) 

and it is moreover easily seen (by setting m = N — 1) that this system is 
consistent with the constraint (Q, provided the parameters aj satisfy the 
following restrictions: 

a 4 = a 6 = 2a 7 — a 12 = a 8 + a 2 = a w = a n = , (26) 
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entailing that the system (I23p becomes 

c m - [(JV - to) a 5 + ax] c m 
+ (JV - 1 - to) (JV + 2 - to) a 7 c m _ 2 - (JV + 1 - to) (JV - 1 - to) a 2 c m _i 

-m [(2 JV - 1 - m) a 9 + a 3 ] c m = ,(27) 

and likewise the system (T2"lj) becomes 

. , , 2a 7 

z n — ax z n -+- a 2 -+- 03 z n 

iV 

~l~ ^ ^ { (^n ^m) [2 "I - (^n "I - ^m) ^5 

m=l 

m^n 

+ 2 (a 7 - a 2 z n + a 9 z*)]} . (28) 

It is easily seen that these systems coincide with ([7]) and after the following 
trivial relabeling of the parameters: aj = Aj for j = 1, 2, 3, 5; a 7 = 
a$ = Aq. And it is moreover plain that the constraint (Q corresponds to the 
constraint ([8]), since the relation ([2]) among the zeros z n (t) and the coefficients 
c m (t) of the polynomial ip clearly entailsj 

N N N 

cN(t) = h n n^w > c *-iw = h" _i e n ' ( 29a ) 



hence 



n=l m=l n=l 



A? 



c^^-c^^-L (29b) 

n=l 

The second possibility is to set a\ 2 = (but an ^ 0) and to then require 
that the system of evolution equations fl23j) be consistent with the constraint 
ffTB"]) . It is indeed plain that the system (T23]) becomes then linear: 

c m - [(N + 1 - m) a 4 - an] c m _i - [(JV - m) a 5 + ai] c m + m a 6 c m+ i 
+ (JV + 1 - to) (JV + 2 - to) a 7 c m _ 2 + (iV + 1 - to) [(JV - to) a 8 + a 2 ] c m _i 

-m [(2 JV - 1 - to) a 9 + a 3 ] c m + m (m + 1) a 10 c m+ i = ,(30) 



2 Actually to reach these conclusions one needs the additional condition cjv(i) ^ 0, 
namely (see (|29at ) z n (t) ^ 0; but the fact that such a condition must be imposed initially, 
and that it shall subsequently hold for all regular solutions of the Newtonian equations of 
motion ([6]), is implied by the presence of the term proportional to in the right-hand 
side of these equations of motion. 
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and it is moreover easily seen (by setting m = N) that this system is consis- 
tent with the constraint (TT5|) . provided the parameters aj satisfy the following 
restrictions: 

CI4 — an = 0*7 — a-2 — (N — 1) ag + 03 = 012 = , (31) 

entailing that the system (123]) becomes 

Cm - (N - m) a 4 c m _i - [(N - m) a 5 + a x ] c m + ma 6 c m+ i 

+ (N + 1 - m) (N - m) a 8 c m _i 
—m (N - m) a g c m + m (m + 1) a w c m+ i = , (32) 

and likewise the system ([21]) becomes 

z n = a x z n - (N - 1) a 9 z n - 2 (N - 1) a 10 z% + a 4 — 

V 

+ { — 2 m) 1 [2 i n Z m + (z n + Z m ) (04 + 05 2 n ) 

m=l 

+a 6 (i„ 2 m + z m z n ) z n + 2 (a 8 z n + a 9 z n 2 + do z„)]} ■ (33) 

It is easily seen that these systems coincide with ffl6]) and ffl~5]) after the 
following trivial relabeling of the parameters: aj = Bj for j = 1,4,5,6; 
a 8 = £> 7 , a 9 = B 2 , aio = -83- And it is moreover plainEl that the constraint 
( TT8]) corresponds to the constraint (TTTj) . see (129aj) . 

3 Behaviors of the TV-body models belonging 
to the first class 

In this section we discuss - in somewhat more detail than done in the intro- 
ductory Section 1 - the behavior of the solutions of iV-body models belonging 
to the first class, see (Q. 

3 Again with the additional condition c/v(i) 0, which need not be highlighted since the 
presence of the term proportional to B4 in the right-hand side of the Newtonian equations 
of motion (|15p guarantees automatically its validity for all regular solutions of this system 
of ODEs. 
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The general solution of the linear system of ODEs with is given 
by the formula 



N 

m = E [ 



(34) 



Here the A^-vectors v^^' are the 2 A/" eigenvectors corresponding to the 2 A^ 
eigenvalues of the (generalized) eigenvalue equation 



(A 2 1 + AM« +M( 2 )) 2 = 



(35) 



where the N x N matrices M} 1 ' and M^ 2 ' are defined componentwise as 
follows (see (J7j)): 

M« = - [(N - m) A, + Ai] , (36a) 



M 



m( 2 ) 

m,m 

(2) 

m,m— 1 



-m [(2JV-m-l) A 6 + A 



3J , 



(iV + 1 - m) {N - 1 - m) A 



2 j 



M, 



(2) 



m,m— 2 



(N — 1—m) (N + 2 - m) A 



i • 



(36b) 
(36c) 
(36d) 



with all other elements vanishing. As for the 2A^ constants they are 
arbitrary except for the two requirements 



N 



< N—l ^ It U N-1 



(37a) 



A/ 







(37b) 



which clearly correspond to the constraints (11 ip via (134)) . Of course, in the 
context of the initial-value problem, the 2 A^ constants are determined 
by the 2N initial data c m (0), c m (0) (including (fill) ) via the system of 2 AT 
linear equations 



N 

,(oi = E[ 

£=1 



(38a) 
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N 



UO) = £ 7i +) ^ )(+) A{ +) + 7 , H «J?W A H . (38b) 



=1 



And in the context of the initial-value problem for the Newtonian iV-body 
model (|SD the initial values c m (0), c m (0) are determined by the initial values 
z n (0), z n (0) via the following relations implied by (J2J): 



V AT 



n=l m=l 
N N N 



- £ z n (0) l[[z- z m (0)} = UO) z N ~ m . (40) 

n=l m=l m=l 

Remark: let us emphasize that, while in the context of the initial-value 
problem for the iV-body problem ([2D the constants 7^ depend on the initial 
data (as just explained), the 2N eigenvalues A^ - characterizing via (JMj) 
the time evolution of the coefficients c m (t) of the polynomial ip(z, t) whose N 
zeros z n (t) yield the coordinates of the N moving particles - do not depend 
on the initial data, but only on the parameters Aj that specify the iV-body 
model ([HD under consideration. 

We now notice that the N x N matrix M} 1 ^ is diagonal, and the N x N 
matrix M} 2 ^ is triangular (Mnm = if n < m). Hence the eigenvalues 
are just the roots of the N (decoupled) second-order equations 

A 2 - [(N-m) A 5 + At) X-m [{2 N-m- 1) A 6 + A 3 ) = , (41) 



1. e. 



A^ = \{(N- m) A, + A x ± A m } , (42a) 

= [(N-m) A 5 + i 1 ] 2 + 4m [{2 N-m- 1) A 6 + A 3 ] 
= (NA 5 + A 1 ) 2 +(Ai-AA 6 ) m 2 

+2 [2A 3 + 2 (2N- 1) A 6 -A 5 (A x + N A 5 )] m . (42b) 

The behavior of the solutions of the iV-body problem ([6]) with (jSJ) is given 
by the time evolution of the N zeros z n (t) of the polynomial ip(z, t), see (J2J), 
whose coefficients evolve exponentially in time as entailed by (|34|) with fT42|) . 
The study of the time evolution of the particle coordinates z n (t) is therefore 
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reduced to the study of the motion of the zeros of a polynomial whose coef- 
ficients depend exponentially on time. In the generic case - characterized by 
exponents not all of which are purely imaginary - the asymptotic behavior 
of these zeros in the remote past and future - characterizing the qualitative 
behavior of the iV-body model under consideration - can therefore be easily 
evinced by the treatment provided in Appendix G (entitled "Asymptotic be- 
havior of the zeros of a polynomial whose coefficients diverge exponentially" ) 
of Ref. [2]. The subclass of iV-body models (jSJ) characterized by parameters 
Aj satisfying the following restrictions, 



[(N-m) f3 + a} 2 -4m [(2JV-m-l) 77 + 7] >0 , m = l,...,N, (43b) 

(where the 4 constants a, (3, 7, rj are all real) - conditions which are necessary 
and sufficient to guarantee that all the 2 N eigenvalues Am ^ are imaginary 
numbers - is instead clearly characterized by the property that all motions 
are confined. The more special subcase in which all the 2 N eigenvalues 
are integer multiples of a single quantity iuo with uo a positive constant, 
uj > 0, so that all the coefficients c m (t) are periodic with the same period 
T = 2 7r / lj, is discussed in Section 5. 

4 Behaviors of the iV-body models belonging 
to the second class 

In this section we discuss - in somewhat more detail than done in the intro- 
ductory Section 1; but rather tersely, to avoid repetitions of developments 
already elaborated in the preceding Section 3 - the behavior of the solutions 
of iV-body models belonging to the second class, see (TT51) . 

The general solution of the linear system of ODEs ffl6l) with ffl8|) is given 
by a formula entirely analogous to (134"|) . and the subsequent developments are 
also analogous to those reported in the preceding section, except that now 
the N x N matrices My' and M} 2 ^ are defined componentwise as follows (see 



Ai = i a , A 3 = 7 , A 5 = i (3 , A e = r] 



(43a) 



flUD): 




m B { 



(44a) 
(44b) 




[(N -m)B 5 + B 1 ] , 



14 



i<Li = -(N-m)B A , (44c) 
M% m+1 =m(m + l)B 3 , (44d) 
M£l = -m(N-m) B 2 , (44e) 



^,Li = (JV + 1 - m) (iV - m) 5 7 , (44f ) 

with all other elements vanishing. As for the 2N constants 7^ , they are 
again arbitrary except now for the single requirement 

tb^^^+VvP^Xt^O, (45) 



which clearly corresponds to the constraint (120]) via (134 

since the N x N matrices M (1) and M (2) are now neither diagonal nor 
triangular, the computations of the eigenvalues cannot be done now in 
explicit form (in the general case, with the coefficients Bj appearing in (1441) 
unrestricted). Therefore a discussion of the actual behavior of the solutions 
of this second solvable model cannot be done here in this general case to the 
same explicit extent as in the case of the first class of models, treated in the 
preceding section. There is however a subcase of this second model for which 
the treatment becomes closely analogous to that of the preceding section. 
Indeed clearly if 

B 3 = B 6 = , (46) 

the two N x N matrices M} 1 ' and M*- 2 - 1 become both triangular, hence the 
corresponding eigenvalue equation becomes again a second-order algebraic 
equation that can be easily solved^. Indeed it is easily seen that if the re- 
striction (j4"6"j) holds, this second-order equation determining the exponents 
Am ^ that characterize the behavior of the second class of models coincides 
with (a subcase of) the second-order equation (T4T]) determining the exponents 
Am^ that characterize the behavior of the first class of models, provided one 
sets 

Ax = B X , A 5 = B 5 , A G = B 2 , A 3 + (N-1) A 6 = , (47) 



4 The N x N matrices M} 1 ^ and M} 2 ^ become triangular also in the case B4 = Bj = 0, 
but whenever B4 vanishes the novelty of the case treated in this paper disappears, so we 
do not pursue this case. 
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hence the discussion given above for the first class of models becomes ap- 
plicable to the second class of models provided they are restricted by the 
condition (j46p . And this includes of course also the more special, entirely 
isochronous, subcase, as discussed in Section 5. 



5 Isochronous versions of the solvable models 

A current definition of isochronous systems (see, for instance, [4]) attributes 
this property to any dynamical system that possesses at least one open (hence 
fully- dimensional) region in its phase space within which all solutions are 
completely periodic (i. e., periodic in all dependent variables) with the same 
fixed period (of course, the period being independent of the initial data, as 
long as they stay within that isochronicity region). The class of isochronous 
systems is vast [4], and it includes a large zoo of systems interpretable as 
A-body problems characterized by autonomous equations of motion of New- 
tonian type. A class of such isochronous systems is characterized by the 
Newtonian equations of motion (see Example 4.1.2-3 in Ref. [4]) 

K 

z = -iuJZ + ^2F { ~ k) (z,z + iujz). (48a) 

k=l 

Here underlined variables indicate A-vectors, u is a nonvanishing real con- 
stant, K is an arbitrary positive integer, and the K (A-vector-valued) func- 
tions F_(~ k ^ (z, z) are required to be analytic (but not necessarily holomorphic) 
in all their 2 A arguments and to satisfy the scaling property 

E_ { ~ k) (cz,~z) = c- k F { - k) (z,z) , k = l,...,K. (48b) 

Indeed it has been shown [4] that these dynamical systems, f j4*8l) . are isochronous, 
possessing an open (hence fully- dimensional) region in their phase space in 
which all their solutions are completely periodic with period 

m 2vr 
T = — 

It is easily seen that the system (j24H belongs to this class, (j48l . with 
K — 1, provided 



(49) 
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d! = -(2N - l)iu, 
a 3 = 2(N - IV 2 , 
a 6 = 0, 

da = -LU 2 , 



«2 = [an - (N - l)a 4 ]ia;, 
a 5 = 2ia>, 
a§ = ia;a4, 
oio = 0, 



(50) 



while the constants 0-4, O7, (In and ai2 remain arbitrary. It is clear that both 
systems, flS} respectively (fT5l) . fall within this class, provided the constants 
Aj respectively Bj featured by these two solvable models satisfy appropriate 
restrictions (whose explicit determination can be left as a simple exercise for 
the diligent reader). 

But here we like to use the more restrictive definition of entire isochronic- 
ity [4] , stating that a dynamical system is entirely isochronous if all its non- 
singular solutions are completely periodic with period T, or possibly with a 
(not arbitrarily large; indeed, generally rather small [12]) integer multiple of 
T - implying that all nonsingular solutions of an entirely isochronous system 
are in fact completely periodic with the same period, which however need not 
be the primitive period for all of them; or, equivalently, that the property of 
isochronicity holds in the entire phase space (with the possible exceptions of 
a lower-dimensional set of initial data yielding singular solutions). 

One can then assert (see Appendix B for a proof) that the iV-body prob- 
lem (ED is indeed entirely isochronous provided the 6 constants Aj it features 
satisfy the following 4 restrictions: 



Ai 
A 3 
A 5 
A. 



(ki - Nk 2 ) 



2 + h) 
k 2 iuo , 



N 



2 ) ( fc 2 - h) 
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k 2 

— u 2 



(51a) 
(51b) 
(51c) 
(51d) 



where u is a positive constant, u> > and the 3 numbers ki,k 2 , k 3 are integers, 
unrestricted (i. e., positive, negative or vanishing) except for the requirement 



fci^-mfc 3 for m = l,...,N . (51e) 

Note that the two constants A 2 and A^ remain completely arbitrary. Also 
note that, if the two integers k\ and k 2 vanish - entailing that A\ and A 5 
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vanish, A\ = A 5 = 0, and A3 = (1 — 2N) A 6 , the only remaining restriction 
on Aq being that it be a positive real number, A 6 > - then the equations 
of motion ([6]) become real (of course provided real values are also assigned 
to the two a priori arbitrary constants A2 and A4). 

But there is a second, different class of entirely isochronous systems, that 
is obtained from (a subclass of) the solvable system §5§ with ([HD via the 
standard "trick" procedure (see for instance Ref. [4]), as we now explain. To 
this end we take as starting point the following special case of the model ([6]): 

/-II 2A 4 ^ CnCm + A/t 

Cn = + 2 2^ 7 _{- ' ( 52 ) 

which corresponds to ([6]) with all the constants Aj set to zero except A4 
(and moreover the purely notational replacement - convenient for what shall 
immediately follow - of the dependent variables z n {t) with the dependent 
variables ( n (r), and of the independent variable t with the independent vari- 
able r, appended primes denoting of course differentiations with respect to 
this new independent variable r). We then set 

/ \ / • \ j- 1 \ ex P (i^*) — 1 /^ oN 
z n (t) = exp (-iut) Qn (r) , r = : , (53) 

1 UJ 

with lo a positive constant, and we thereby obtain the new iV-body model 
characterized by the {autonomous) Newtonian equations of motion 

2A 4 

z n — — 1 UJ z n 

N 

+2 ^2 {(z n - z m )~ l [(z n + \uz n )(z m + \uz m ) + A 4 }} . (54) 

771=1 
771^77 

Due to the way this model has been obtained it is clear that it describes 
an entirely isochronous iV-body problem - provided its initial data are con- 
strained by the conditions (TlOl) . guaranteeing the validity of (IE]) (both for this 
model and for the model fl52l) ). It is also easy to verify that this model is not 
a special case of the entirely isochronous iV-body problem identified above, 
i. e. §6§ with (I5T1) (and of course with ©). 

Let us turn now to a discussion of entirely isochronous variants of the 
second class of solvable models treated in this paper. The identification of 
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such models cannot be explicitly achieved in the general case along the same 
lines as done above - in the initial part of this section - for models belonging 
to the first class of solvable systems^ because in this second case one cannot 
obtain explicitly the exponents Am , as explained in the preceding section 
(see the paragraph after (135]) unless the restriction (j46p holds, in which case 
the treatment given in the first part of this section for the first class of system 
can be immediately extended to the second class via the relations (HT|) . which 
however entail (see the last of these equations (|47|) ) that the restrictions (I5T1) 
must now be complemented by the additional restriction 

k 2 + k 3 = or 2 k x = N(k 2 - k 3 ). (55) 

It is moreover again possible, also for the models of the second class, 
to proceed via the standard trick procedure. Since this approach is quite 
analogous to that described immediately above for the first class, we deal 
with it quite tersely. 

The starting point is now the solvable system ffl5l) with all constants 
vanishing except B4, which we now write as follows: 

_ £^ C _|_ 2 C Cm + -^4 (Cn + Cm) ^g-j 

Cn 1 , Cn Cm 

m=l,m^=n 

It is now again appropriate to use the version (|53|) of the trick, obtaining 
thereby the iV-body problem - clearly entirely isochronous, because of the 
way it is obtained - characterized by the following equations of motion: 



. . z n + iuz n 

Z n = —lLUZ n + -D4 

+ 2 fat + iuzn) (An + iuz m ) + B± [z n + z m + iu (z n + z m )} 

m=l Zn ~ Zm 



6 Numerical results 

In this Section we perform a numerical integration of the equations of motion 
([6]) of the first model in order to illustrate the previous findings. The nu- 
merical integration has been done with an embedded Runge-Kutta method of 
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order 8(5,3) with automatic step size control, as developed by Prince and Dor- 
mand [13]. The integration and the graphical output have been performed 
with the software Dynamics solverO developed by J. Aguirregabiria. 

In order to simplify our numerical treatment we have set in ([6]) N = 3, 
and A\ = A 2 = A 5 = 0. Hence, our numerical study focuses on the following 
system of three nonlinear complex ODEs of second order: 

z n = A 3 z n h - Zm)' 1 [2 z n z m + 2(A A + A 6 z%)]. (58) 

*M -. 

m=l 

Generic (chaotic) behaviour 

When A4 = this system is solvable and it is a particular case of the family 
discussed in Chapter 2.3.3 of [2]. When A4 7^ and no restrictions are 
imposed on the initial data, the system (1581) corresponds to the motion of the 
zeros of a 3 rd degree polynomial whose coefficients evolve in time according 
to the nonlinear system fl23|) . The system behaves in general in a chaotic 
manner as it can be evinced from the numerical results in Figures le, If and 
2. The initial data for that integration have been chosen to be 

Zl (0) = 1.5 - 2i, ii(0) = 2-i 

z 2 (0) = -l + 0i, i 2 (0) = -l + i (59) 
z 3 (0) = 1.2 + 0.4L £3(0) = 2.32 - 0. 76i 

which do not satisfy the constraints ffTUl) . while the values of the coefficients 
are 

A 3 = -5.1tt 2 , A 4 = 15, A 6 = n 2 . (60) 
Quasi-periodic behaviour 

However, for the same values of these parameters, flBDl) . initial data can be 
chosen within the algebraic submanifold defined by ffTUl) . One such assign- 
ment is 

*i(0) = 2-2i, ii(0)=2-i, 

z 2 (0) = -l + 0i, i 2 (0) = -l + i, (61) 

23(0) = 1.2 + 0.4i, i 3 (0) = 2.32 - 0.76L 

5 This program is available at http://tp.lc.ehu.es/jma/ds/ds.html 
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As discussed in Section 2, now the system (I58|) is solvable, and it corresponds 
to the motion of the zeros of a monic 3 rd degree polynomial whose coefficients 
evolve in time according to the following linear system: 

c\ - (4A 6 + A 3 ) Cl = 0, (62a) 
c 2 = 0, (62b) 
c 3 - 2A A c t - 3(2A 6 + A 3 )c 3 = 0. (62c) 

The solution of this system is trivial and from (1431) and the discussion in 
Section 4, we see that it suffices that 

AA Q + A 3 < 0, 2A 6 + A 3 < 

for the orbits to be confined. As it follows from (1541) with (1421) . the coefficients 
of the polynomial are in general quasi-periodic functions of time, hence also 
the behaviour of the system (TTOT) is quasi-periodic as it can be seen in Figures 
lc, Id and 2, where we have again assigned the values (IBTJI) of the parameters, 
and the same initial data (IBTj) . 

Isochronous behaviour 

Following the discussion of Section 5 and Appendix B, the family of solvable 
problems discussed above includes a subclass of entirely isochronous systems. 
They occur when all the frequencies of the system are commensurable. A 
look at (joTj) reveals that these conditions are met if we choose 

A 1 = A 5 = 0, A 3 = -- A klu\ A 6 = ^k 2 3 u 2 . (63) 

We have set k 3 = 1 and uj = 2ir so that the fundamental period is T = 1. 
The choice of parameters is then 

A 3 = -5tt 2 , A 4 = 15, A 6 = ix 2 . (64) 

and the initial data are the same as in the previous case ( I6TT) . The resulting 
motions are shown in Figures la, lb and 2. 

7 Outlook 

The novel cases treated in this paper do not exhaust all the possible solvable 
many-body problems associated to the motion of the roots of polynomials 



21 



whose coefficients satisfy a linear system of ODEs. This question is related 
to the classification of exceptional subspaces with co-dimension k > 1, 
which has not yet been completed. However, as mentioned in Appendix A, 
some results based on X 2 spaces and operators will be reported soon [5]. 
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Appendix A 

In this Appendix we explain the general method to write down certain PDEs 
in 1+1 dimensions possessing solutions which are monic polynomials ip{z,t) 
of degree n in z whose N coefficients are time-dependent functions c m (t). 
The requirement to ensure the solvability of the nonlinear iV-body problem 
characterizing the time evolution of the N zeros z n (t) of such a polynomial 
is that the coefficients c m (t) evolve according to a linear system of ODEs. 

Let us denote by Vn the vector space of polynomials in z of degree less 
than or equal to N: 

V N = span {l,z, z 2 ,...,z N }, dimV N = N + 1. (65) 

We shall denote also by T>2(Vn) the vector space of differential operators of 
order 2 in z with analytic coefficients that leave the space Vn invariant, i.e. 

V 2 (V N ) = {T = a 2 (z)D zz + ai (z)D z + a {z) | TV N C V N } . (66) 

Notation: here and below we use the notation D z = d/dz, D zz = d 2 /dz 2 , 

D t = d/dt, D tt EE d 2 /dt 2 , D zt EE d 2 /dzdt. 
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Consider the action of the most general linear second-order partial differ- 
ential operator in z and t: 

L = aiD a + a 2 D zt + a 3 D zz + a 4 D z + a 5 D t + a 6 , a { = a* (2, t) (67) 

on the following time-dependent polynomial of degree N in z (see fl2])), 

N N 

lf>(z,t) = 1[[Z- Zn(t)) =Z N + Ys C m(t) Z"'™ ■ (68) 
n=l m=l 

The restriction to second-order operators ensures that the nonlinear many- 
body problem for the zeros will have at most two-body interactions, which is 
the more physically interesting case. If many-body interactions are allowed, it 
is obvious that the class of such solvable systems that can be constructed with 
this method is larger. Since we want the evolution of the coefficients c m (t) to 
be a linear autonomous (hence solvable) system, we assume a,i(z,t) = ai(z) 
and without loss of generality we also set 01(2) = 1. 

We now examine the conditions on a^z) such that the equation L[ip] =0 
implies a linear system for c m (t). These conditions arise from imposing that 
L acting on i/;(z, t) produces a polynomial of degree N — 1. They are satisfied 
by a second- order differential operator in z and t only in the following cases: 

• the second-order differential operator D a ] 

• one of the following 4 differential operators of first-order in t: 

D t , D zt , zD zt , z[zD zt -{N-l)D t ] ; (69) 

• a second-order differential operator in z that maps Vn to Vn-i- 

The characterization of vector spaces of linear (and nonlinear) differential 
operators in one and several variables of any given order that map Vn to 
Vn-ic has been treated in [9] . Working out all the possible cases for a second 
order differential operator in one variable that maps Vn to Vn-i is a simple 



exercise that produces the following result: 

T+ 1 = z(zD z -N+ l)(zD z -N + 2), (70a) 

T 2 ° = z 2 D zz - N(N - 1), (70b) 

17 1 = zD zz , (70c) 

T 2 2 = D zz , (70d) 

If = zD z - N, (70e) 

Tr 1 = D z . (70f) 
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A linear combination of the eleven operators written above is precisely equa- 
tion ([TJ) of this paper, which was treated in Section 2.3.3 of [2] (and several 
specific cases were then investigated in Section 2.3.4). However, the question 
arises whether this is the most general class of PDEs in z and t which ad- 
mit as solutions polynomials in z whose coefficients are functions of t that 
evolve according to a linear system. The recent discovery of the so-called 
exceptional polynomial subspaces of Vn shows that other PDEs with those 
properties exist (and therefore other solvable many-body problems) provided 
some constraints are imposed. 

An exceptional polynomial subspace Ai^ of co-dimension k in Vn is 
defined by the property that some second-order differential operators that 
preserve M$ do not preserve Vn- More specifically, consider a space M$ C 
Vn generated by N + 1 — k linearly independent polynomials, all of them of 
degree at most iV: 

M ( n = span{pi(z), . . . ,p N +i-k{z)} ■ (71) 

We will say that M N * s an exceptional polynomial subspace of co-dimension 
k in Vn (for short, an X^-space) if T>2(A4$) ^(Vn)- Such exceptional 
subspaces do exist, and they provide novel differential operators with poly- 
nomial eigenf unctions. Exceptional subspaces have been first analyzed in 
the context of quasi-exactly solvable potentials in quantum mechanics [7,10]; 
they are also connected with the Darboux transformation [8] and with non 
classical families of orthogonal polynomials [11]. Xi-spaces have been fully 
classified arriving at the result that there is essentially one such subspace up 
to projective transformations. This space is precisely 

X 1 = sptm{l,z 2 ,z 3 ,...,z N } , (72a) 

which can also be characterized as 

X 1 = {peV N \p'(0) = 0} . (72b) 

It can be shown that dimV 2 (X 1 ) = 7. The most interesting element of this 
space is the operator 

D zz - 2 -D z , (73) 

which preserves X\ but not Vn- The first class of models treated in this 
paper - see (jSJ)- is related to this exceptional space X±, and the novel term 
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proportional to the constant A4 in the many-body problem ([H]) is precisely 
that due to the inclusion of this operator. The constraint on the many-body 
problem (jSj) is precisely the one that defines the exceptional subspace X\. 
The second class of many-body models treated above has a different origin, 
not related to the exceptional subspaces of Pat, being instead associated with 
the role of the time-derivative, see (161)1) . 

Exceptional subspaces of higher co-dimension exist but a full classifica- 
tion is not yet available. Some new many-body problems with constraints 
associated to X 2 -spaces will be treated in a forthcoming publication [5] , their 
main novelty being that the coefficients of the new differential operators are 
not only inverse powers but rational functions of z. 



Appendix B 

In this Appendix we justify the assertion made in Section 5, that the condi- 
tions (l5"Tj) guarantee that the iV-body model ([6]) with (1101) is entirely isochronous. 

Indeed, clearly this iV-body problem is entirely isochronous if all the 
solutions of the corresponding linear problem ([7j) are completely periodic with 
the same period T (readers for whom this is not clear are advised to consult, 
say, Ref. [4]). Clearly a necessary and sufficient condition for this to happen is 
that the 2N eigenvalues Xm \ see (l4"2"j) . be all integer multiples of a common 
imaginary constant, say 

)$ = kghu>, (74a) 

with uj a positive constant, u = 2n/T > 0, and the coefficients km all 
integers, and moreover that 

AL +) + A« , (74b) 

for all values of the index m, i. e. for m = 1,...,N. Clearly the first of 
these two relations, (I74al) . can only be true if the quantity A^, see (142bl) . is 
a perfect square for all values of m, so that 

A m = a + (3m, (75) 

clearly entailing (see f!42bj) ) 

a = NA s + A 1 , (3 2 = A^-4A 6 , (76a) 
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and (as we hereafter assume) 

2A 3 + 2 (2N- 1) A e - A 5 (N A 5 + Ax) = a f3 . 
Via f)42p these formulas yield 

Ai +) = « + | (J3-Az) , AH = -|(/? + A) 

Hence f )74al) entails 

a = i a; , (3 — A 5 = 2 j iu , j3 + A 5 = 2kiu 



(76b) 



(77) 



(78) 



with ki,j,k integers, and the (sum and difference of the) last two formulas 
entail 

A 5 = k 2 iu, f3 = k 3 iu, (79) 
with k 2 and k 3 integers. Hence the two equations (j 76 aft yield 



Ax = (k, - Nk 2 ) iu , A 6 = fc ' , k * u 2 



and from (176bj) one finally gets 



A 3=2 ^ 2 + k ^ 



4 



N - - ) (k 2 - k 3 ) - h 



(80) 



(81) 



The 4 formulas (T5TT) are thereby proven. And it is moreover clear that the 
condition (I74bj) entails the requirement (15 lei) , thereby completing the proof 
of the results reported in Section 5. 
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Figure 1: Left column: trajectories of the 3-body system ( 158|) in the pe- 
riodic (la), quasi-periodic (lc) or chaotic regimes (le). Right column: 
Stereographic projection of the trajectories: the positions of the three parti- 
cles are plotted at every integer multiple of T = 1 up to 10 3 iterations. 
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Figure 2: A plot of Ini2 2 (t) for the choice of parameters and initial data 
described in Section 6. The periodic, quasi-periodic and chaotic behaviours 
are manifest. 
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